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We calculate the dynamical conductivity of AA-stacked bilayer graphene as a function of frequency 
and in the presence of a finite chemical potential due to charging. Unlike the monolayer, we find 
a Drude absorption at charge neutrality in addition to an interband absorption with onset of twice 
the interlayer hopping energy. At finite doping, the interband absorption exhibits two edges which 
depend on both chemical potential and interlayer hopping energy. We study the behaviour as a 
function of varying chemical potential relative to the interlayer hopping energy scale and compute 
the partial optical sum. The results are contrasted with the previously published case of AB- 
stacking. While we focus on in-plane conductivity, we also provide the perpendicular conductivity 
for both AB and AA stacking. We also examine conductivity for other variations with AA-stacking, 
such as AAA-stacked trilayer. Based on proposed models for topological insulators discussed in the 
literature, we also consider the effect of spin orbit coupling on the optical properties of an AA-stacked 
bilayer which illustrates the effect of an energy gap opening at points in the band structure. 
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I. INTRODUCTION 

Graphene continues to provide a rich platform for in- 
vestigations into the physics of massless Dirac fermions. 
Initially studies of graphene were limited to the realm 
of theory where the low energy linear dispersion^ and 
chiral nature of the honeycomb carbon lattice were 
shown^ to result from a simple nearest-neighbor-hopping 
tight-binding Hamiltonian which at low energy maps 
on to a Dirac Hamiltonian for massless fermions with 
Fermi velocity vf- With the experimental realization 
of graphene^, a considerable literature has now accu- 
mulated which has uncovered a variety of exotic effects, 
such as an unusual quantum Hall effecl£i£, giant Faraday 
rotation 7 , plasmarons 8 , and so on, some of which has 
been summarized in reviews^—. 

Bilayer graphene is also of intense interest as it too 
shows an unusual quantum Hall effec t 13 ' 14 and indeed its 
low energy tight-binding Hamiltonian maps to an equa- 
tion for chiral fermions with an effective mass^ based on 
an interlayer hopping parameter 7. In addition, it has 
been seen that bilayer graphene can develop a sizeable 
band gap which is tunable by charge doping J£ Recent in- 
terest in bilayer graphene physics has focused on the large 
degeneracy at the charge neutrality point which pro- 
vides opportunity for instabilities leading to new ground 
states. See Ref for a summary of the literature on this 
point and also a general review of the properties of bi- 
layer graphene. The natural form for bilayer graphene 
is the so-called Bernal or AB-stacking which is the ba- 
sis of the parent compound graphite from which it is 
usually derived. Consequently, past work has primarily 
focused on this stacking configuration. However, more 
recently, Moire patterns seen in scanning tunneling mi- 
croscopy imaging of graphene bilayers and multilayers 
point to alternative stackings where one layer is rotated 



by some angle relative to the otherj 17 i 18 This is sometimes 
referred to as twisted or misaligned bilayer graphene. In 
these systems, the electronic properties are modified at 
low energy such that monolayer behaviour appears along 
with a reduced Fermi velocity^ In these systems, it is 
possible to have regions which are rich in AB-stacking 
and regions which display mainly AA-stacking. These 
types of stackings are shown in Fig. [Ha) and (b) , respec- 
tively. Here, A and B refer to atoms on the two triangu- 
lar sublattices of the honeycomb lattice and the stacking 
is in reference to whether the A atom of one plane is 
stacked over the A or B atom of the other plane. For 
Bernal AB-stacking only half the atoms are aligned on 
top of each other and the other half sit over the center 
of the hexagon in the opposite layer. For AA-stacking 
all atoms are matched up between the two layers. For 
very small twist angles, the regions of AA-stacking have 
been suggested to provide localization effects Very 
recently, AA-stacked graphene has also attracted inter- 
est due to research which has identified such stacking in 
certain samples, potentially making this another experi- 
mentally accessible system to study 21 ' 22 For AA-stacking 
there is also the prediction for new ground states to oc- 
cur, such as antiferromagnetism 2 ^. 

As a result, we are motivated by these develop- 
ments to examine the dynamical conductivity of doped 
AA-stacked graphene to elucidate features which would 
demonstrate unique properties of this system and allow 
for the identification of characteristic energy scales asso- 
ciated with the band structure. Moreover, as the optical 
properties of graphene are of considerable importance for 
technological applications, all variants of graphene are 
also of potential interest and should be examined. The 
dynamical conductivity of graphene has been extensively 
studied theoretically^— and experiments have largely 
verified the expected behaviour— ~—. Likewise, the con- 
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FIG. 1. (Color online) (a) Two AB-stacked graphene sheets 
with the dark (blue) dots representing one sublattice and the 
light (green) dots, the other, (b) Two AA-stacked graphene 
sheets. 



ductivity for Bernal-stacked bilayer graphene has been 
predicted theoretically^— and observed^— . There 
has also been work on magneto-optical conductivity of 
graphene in which theory and experiment are also in good 
agreement. Indeed, a review of this literature may be 
found in Ref4i. Some preliminary work on the absorp- 
tion coefficient of undoped AA-stacked graphene in zero 
magnetic field has been reported^ however most materi- 
als naturally occur with charge doping where the Fermi 
level or chemical potential fi is away from charge neu- 
trality (// = 0). Furthermore, the interesting feature for 
practical applications is the variation of optical proper- 
ties with doping, usually achieved through a field effect 
transistor structure. 

In the following, we provide a thorough examination 
of the finite frequency conductivity for AA-stacking, for 
both in-plane and out-of-plane response. Unlike a single 
monolayer, AA-stacked graphene shows a Drude response 
in the in-plane conductivity at charge neutrality along 
with Pauli blocking at low frequencies below the onset 
of a flat interband absorption. This interband absorp- 
tion splits at finite doping into two interband absorption 
edges leading to flat universal values associated with one 
and two layers. In terms of the interlayer hopping en- 
ergy 7, the partial optical sum, which is a measure of the 
transfer of spectral weight, shows distinct behaviour for 
/i < 7 versus jj, > 7. The perpendicular conductivity has 
a strong response at 27 at all dopings. This is contrasted 
with the case for AB-stacking which we also show here as 
we can also provide an analytical formula for this quan- 
tity to add to the literature. Indeed, we have provided 



analytical formulae in almost all cases and physical un- 
derstanding of our results are given. We also contrast 
the AA-stacked case with that for AAA-stacking. Fi- 
nally, because of the connection between the Dirac na- 
ture of graphene and topological insulators (TIs), we also 
follow-up on a toy-model^ by providing the conductivity 
for two AA-stacked sheets with spin orbit coupling in one 
or both planes, potentially mimicking weakly coupled TIs 
or a TI in proximity to a metallic sheet. 

Our paper is organized as follows. In Section II, 
we review our theoretical calculation for the dynami- 
cal conductivity of AA-stacked graphene. Our presen- 
tation follows that for AB-stacked bilayer graphene done 
by Nicol and Carbotte^ which uses many-body Green's 
function which easily allows for further theoretical de- 
velopment, such as the inclusion of a self-energy from 
impurities^, electron-phonon interaction^—, electron- 
electron interactions^— , etc. In Section III, we dis- 
cuss the results of the AA-stacked case, examining both 
in-plane and perpendicular conductivity and contrasting 
with the AB-stacked case. We discuss the effect of biasing 
the bilayer. We also consider the theory for other varia- 
tions on the AA-stacked case in the subsequent sections. 
For instance, we examine the case for the AAA-stacked 
trilayer in Section IV and report results for models with 
spin orbit coupling in AA-stacked bilayer in Section V. 
Our conclusions arc found in Section VI. 



II. THEORY FOR AA-STACKED BILAYER 

To derive the optical conductivity of AA-stacked bi- 
layer graphene, we follow the method shown in the work 
of Nicol and Carbotte^ for the case of AB-stacked bi- 
layer graphene. This is based on the Kubo formula for 
the current-current response function and the many-body 
Green's function approach^ Thus, to begin we must first 
examine the band structure and provide an expression for 
the electronic Green's function. For the case of AA stack- 
ing, an A (B) atom in the upper layer is stacked directly 
above an A (B) atom in the lower layer, see Fig. [TJb), 
as opposed to the typical Bernal stacking shown in Fig. 

ma). 

For AA stacking, the single spin Hamiltonian is given 

by 

H = -*£ + H.c.) -tJ2 (f 2n + s a 2n + H.c.) 

rfc,5 n,S 

+ lY,{ a 2n a lr l +b\ n b lrl +H.C^ . (I) 
n 

The first two terms are the nearest-neighbour intralayer 
hopping terms for electrons to move within a given plane 
with hopping energy t ~ 3 eV. The two planes are in- 
dexed 1 and 2. As a consequence of the geometry of 
the honeycomb lattice, each sheet has two inequivalent 
atoms labelled A and B. The operator ai n annihilates 
an electron which is on an A-atom site with site label n 
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in the graphene sheet indexed by i. The label n indexes 
the sites of the triangular Bravais lattice. Conversely, 
b\ n+s creates an electron in sheet i on the neighboring 
site at the position n + 6, where S is one of three possible 
nearest-neighbour vectors given by Si = — (ai + a2)/3, 
<$2 = (2ai— a2)/3 and S3 = — (ai — 2ei2)/3. The primitive 
vectors of the triangular sublattice are a\ = (aV3/2,a/2) 
and a 2 = (av / 3/2, — a/2), where |ai| = |a 2 | = \/3a cc 
with a cc the shortest carbon-carbon distance. The third 
term in Eqn. ([1} corresponds to the interlayer hopping 
between graphene sheets. The hopping parameter be- 
tween an A (B) site in one layer and the nearest A (B) 
site in the other layer is given by 7 which is reported to 
be about 0.2 e V 42 ' 54 , which differs in AB-stacked bilayer 
graphene where it is closer to 0.4 eV. There is also a pos- 
sibility to hop between an A (B) site in one layer to a B 
(A) site in the other layer; however, these hopping ener- 
gies are very smallS 4 ^ and thus ignored in our model. 
The Hamiltonian given in Eqn. ([T]) transforms to k space 
in the usual way^ and can be written in the following 
matrix representation: 



H 



7 f(k) 

f*(k) 7 

7 f(k) 

/*(fc) 7 



(2) 



where f{k) = -tJ2s 



and we have used the eigen- 



vector \l/ = (ai k, i>2ki a<2fc, bi k) following the notation of 
McCann 56 . The band structure is given by the eigen- 
values of this matrix. Reflecting the fact that there are 
now four atoms per unit cell, we obtain the following four 
energy bands: 



e a (k) = ±[\f(k)\ + (-r 7 ], 



(3) 



where a — 1 and 2 and |/(fe)| is the energy dispersion for 
a single sheet of graphene. We essentially have two copies 
of the band structure of monolayer graphene, one shifted 
by —7 and the other by +7, or bonding and antibonding 
bands, and indeed we will see that this provides part 
of the physics which enters the dynamical conductivity. 
As our interest is to understand the conductivity at low 
energies, we choose to expand f(k) around the K point 
of the Brillouin Zone to obtain f(k) — hvpke 1,8 , where 
vp = \/3ta/2h and 6 is the fc-space angle around the 
K point. The low energy band structure can be seen in 
Fig.[5]where it is compared to that of the familiar Bernal 
stacking. As the physics of the conductivity associated 
with the K' will be the same as for the K point, it is 
sufficient to work only about the one K point in what 
follows and multiply the result by a factor of two for the 
so-called valley degeneracy associated with the two K 
points per unit cell. 

We can also provide, as others have shown^i, an ana- 
lytic expression for the total double spin density of states, 
N(e): 
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FIG. 2. (Color online) Top: Low energy dispersion for a bi- 
layer with Bernal AB-stacking (left) and AA-stacking (right). 
Bottom: Low energy density of states [in units of 2y /tt(Tivf) 2 ] 
for AB- and AA-stacked bilayers, left and right frames respec- 
tively. 



which results from the sum of two Dirac cone density of 
states shifted relative to each other by 27. A plot of the 
low energy density of states in units of 2j/n(hvF) 2 for 
AA-stacked bilayer graphene is shown in Fig. [2] and is 
contrasted with that for AB-stacking^. 

With our Hamiltonian, it is straightforward to deter- 
mine the Green's function G(z) from G~ l {z) — zl — H . 
Thus 
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The only elements of the Green's function that contribute 
to our final expressions for longitudinal and perpendicu- 
lar optical conductivity are Gu, G12, G13 and G14. We 
will only show these elements explicitly: 



Gxx{z) 



G 12 (z) 



G 13 (z) 



z 3 ~z( 7 2 + |/(fc)| 2 ) 
[z - e 1 )(z + £i)(z - e 2 ){z + e 2 ) ' 

2 7 */(fc) 

(z - Ei)(z + £i)(z - £ 2 ){z + £ 2 ) ' 



z 2 7 + 7 |/(fc)| 2 -7 3 

[Z - £l){z + £i)(z - £ 2 ){z + £2) ' 



(6) 



(7) 



(8) 



4 



and 

r ( x z 2 f(k) +1 2 f(k)-f(k)\f(k)\ 2 
14[Z) ~{z- ei)(z + e 1 ){z - e 2 )(z + e 2 ) ' W 

The finite frequency conductivity is calculated through 
the standard procedure of using the Kubo formula^, 
where the conductivity is written in terms of the retarded 
current-current correlation function. From this the real 
part of the conductivity can be written as2£ 



2Q 



^[/(u;- M )-/(u, + Q- M )] 



d 2 k 
{2-Kf 



Tr 



v a A{u + Sl,k)vpA(w,k) , (10) 



where we have used the spectral function representation 
of the Green's function, 



Gij{z)= — 



dto Aij(uj) 



2ir z- 



(11) 



Here a and j3 represent the spatial coordinates x,y,z, Nf 
is a degeneracy factor, f(x) = l/[exp(x/T) + 1] is the 
Fermi function for temperature T and fi is the chemical 
potential taken to be positive here but to accommodate 
for negative values, // just needs to replaced by ev- 
erywhere. Note that we will usually take % = 1 when re- 
ferring to the relationship between energy and frequency 
and restore it when necessary. For our results, we show 
only the T = case. For the longitudinal in-plane con- 
ductivity, cr xx (£l), v a — v/3 = v x where 



v F 

v F 

v F 

v F 



(12) 



The velocity operator can be evaluated from a Pcicrls 
substitution as demonstrated in Ref^ or from hv x — 
dH/dk x . We can then evaluate the trace, drop the terms 
that will vanish upon averaging over angle and obtain an 
expression dependent on the two spectral functions An 
and A13. In the zero temperature limit, the real part of 
the longitudinal conductivity a xx (fl) is then 



<J XX {®) 



2n 



du) 
2^ 



d 2 k 



Avi 



Au(u + n)An(uj) + A 13 (uj + n)A 13 (uj) 



(13) 

In keeping with the low energy expansion of f(k) about 
a single K point, the integral over k, which in general 
is over the first Brillouin Zone, is now taken as an inte- 
gral over a single K point. The degeneracy factor is thus 
Nf = g s g v , where g s = 2 to account for the sum over 
spin which has been ignored up until now and g v — 2 to 



account for a sum over the K and K 1 points of the Bril- 
louin Zone. Furthermore, the upper limit of the k integral 
is taken to be a large cutoff value typical of momentum 
associated with the large bandwidth. It is convenient to 
scale Eqn. (IT3l) by the constant background conductiv- 
ity of a single sheet of graphene^ given by <ro = e 2 /Ah. 
All that remains before we can calculate our conductiv- 
ity is to specify the necessary spectral function elements. 
Given our expressions for G\\ and G13 and Eqn. (jlip . we 
obtain 



An 

and 
A 13 



[6(w - £1) + 8(u + ex) + 5{lu - e 2 ) + 5(uj + e 2 )] 

(14) 



7T 

- [6(cj + £1) - 6(u) - £1) + 5{uj - £ 2 ) 



S(lJ + £ 2 )}. 

(15) 

For our numerical work, we use the Lorentzian represen- 
tation of the delta function, S(x) = (rj j V) / '[rj 2 + x 2 ], with 
a broadening of 77 = O.OI7. The broadening is manifest 
in the optical conductivity as an effective transport scat- 
tering rate of 1/ Ti mp — 2ry due to the convolution of the 
two Lorentzian functions in the conductivity formula. 

We can also examine the perpendicular conductivity, 
a zz (Q), associated with transport perpendicular to the 
graphene sheets. In Eqn. (fTU|). our velocity operator is 
now (v a =vp =v z ) 



( v* 

v 

v 

V v* 



(16) 



where v = ijd/h with d the interlayer distance, d 
is about 3.6 A and 3.3 A for AA- and AB-stacking 
respectivel y 42 ' 54 . This leads to the real part of the zero 
temperature perpendicular conductivity: 



N f e 2 
2VL 



fi — Q 



duj 
2^ 



d 2 k 

(27T) 2 



4H 2 



A u {lu + n)A n (u) ~ A* n {w + Sl)A 12 {u) 



- A 13 (uj + Q)A 13 (oj) + Al^cu + fi)Au(a>) 

(17) 

with An and A i3 given by Eqns. (|T4^) and (fT5|) . respec- 
tively, and 



A12 = 



2s 



[—5(uj - £1) - 5(uj + £1) + S(u> - £ 2 ) 
+ 5(uj + £ 2 )] 



(18) 



and 



Au = - [S{u - £1) - 8(oj + £1) + 5(u> - £ 2 ) 



2£ 

With£=|/(fe)| 



5(u + e 2 )], (19) 
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FIG. 3. (Color online) (a) Real part of the longitudi- 
nal frequency- dependent optical conductivity a xx (Q) of AA- 
stacked bilayer graphene normalized to that of a single 
graphene sheet for the case of /i < 7. Onset of absorption 
at 2|7 — fi\ and 2(7 + /x) correspond the availability of allowed 
transitions within the AA band structure shown in the in- 
set. Inset: Low energy band structure around the K point 
for the case of fi < 7. Transitions can only occur within the 
same colored bands. For the blue bands, the lowest transition 
between the upper and lower cones occurs at Q, — 2(7 + fj,), 
shown by the blue arrow, and for the red bands, the lowest 
transition between cones occurs at Q — 2(7-/1), shown by 
the red arrow. Below these frequencies, only intraband tran- 
sitions within the same cone for states about the chemical 
potential occur, (b) The case of /j, > 7. 



III. RESULTS FOR AA-STACKED BILAYER 

Here we present results for the longitudinal conductiv- 
ity which is obtained by evaluating Eqn. (|13[) numerically 
using the Lorentzian form in place of the delta function 
and taking the broadening parameter 7/ = O.OI7. Fig. [3] 
shows curves for the case of fj, < 7 and fx > 7 (top 
and bottom frames, respectively). For the case of charge 
neutrality (/j = 0), the conductivity displays a Drude 
response at low frequency due to intraband transitions 
and a flat interband absorption which commences at 27. 
This is quite unlike the case of monolayer graphene which 



would have had a flat interband response at all frequen- 
cies. Here the response is not simply that of two mono- 
layers as one might naively think. This is because the 
AA-stacked bands are essentially two decoupled graphene 
monolayer bands which represent bonding and antibond- 
ing bands and are shifted relative to each other— (as em- 
phasized in the inset where the bands are identified with 
different colors). The matrix elements for the longitudi- 
nal conductivity only allow for transitions between like- 
colored bands.— The transitions must be vertical as the 
photon is a momentum q ~ probe. At charge neutral- 
ity, the minimum interband transition that is not Pauli 
blocked is then 27 and intraband transitions can also now 
occur due to the Fermi level being located away from the 
Dirac point of the monolayer band, this latter feature is 
not present in previous work 42 . This emphasizes that 
knowledge of the band structure as shown in Fig. [2] is 
insufficient (or may be misleading) to the determination 
of the allowed absorption transitions and that the ma- 
trix elements which know about effects of chirality and 
bonding/antibonding are important as well. 

For finite chemical potential, the Drude persists but 
now there are two pieces to the interband absorption 
which onset at 2|7 — fi\ and 2(7 + /i). The behavior is 
different for fi < 7 versus \i > 7. For the case of pb < 7, 
the Drude conductivity remains completely unchanged 
and its weight is set by the value of 7. The interband 
edge that was at 27 in the /i = case is now split into 
two edges moving to lower and higher frequency and as- 
sociated with the onset of allowed transitions for a single 
monolayer band structure shifted up or down by 7, re- 
spectively. This is emphasized by the AA-stacked band 
structure shown in the inset of Fig. EJa) where transi- 
tions can only occur within the same colored bands. For 
the monolayer dispersion shifted down by 7 (blue bands) 
the lowest interband transition from an occupied state 
to an unoccupied state occurs at f2 = 2(7 + /x). For the 
monolayer dispersion shifted up by 7 (red bands), the 
lowest transition between cones occurs at ft — 2\j — /zj. 
Thus at low enough energy, the finite frequency conduc- 
tivity displays the universal background absorption of a 
monolayer <jq but at higher frequency there is a step-up 
to a flat universal background at 2<ro and the transition 
between these steps is tunable with the charge doping. 

For the case of [i — 7 as shown in Fig. EJb), the lower 
edge has disappeared and there is a background conduc- 
tivity of (To for < 47 and 2ctq for Q > A-y. However, in 
this case, the Drude component still remains at very low 
frequency. For /1 > 7 as shown in Fig. EJb) , the lower 
edge reappears, showing the double step in universal con- 
ductivity value and now both edges move to higher fre- 
quency with increased /i. As the area of the conductivity 
is conserved, the lost weight at finite frequency is found 
in the Drude which now increases with /i as one finds in 
the monolayer case. These characteristic features of the 
AA-stacked bilayer are quite different from the case of 
Bernal stacking^ and are not at all the expectation of 
twice the monolayer conductivity either. The presence 
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of the Drude at charge neutrality in the AA-stacked case 
is different from the monolayer and AB-stacked bilayer 
where no such feature exists for /i = 0. These special 
features of AA-stacked graphene might prove useful for 
applications where optical response is tuned by doping 
(or gating) to be like a switch with three settings: off or 
0, on at half setting (<7o) and on at full setting (2<7o). Use 
of tuning by gating has been demonstrated for graphene 
terahertz modulators where the intraband transitions are 
used in this case.— Another significance of the result is 
that for the spectral range between 2|/i — 7I and 2(7 + 
one might not be able to separate monolayer from AA- 
stacked bilayers by optics alone. Overall, the dynamical 
conductivity is quite distinct from that of AB-stacked bi- 
layer where no such steps occur and the conductivity is 
not flat in the low frequency spectral rangei 34 i 35 

These results are embodied by a closed algebraic for- 
mula for the real part of the longitudinal conductivity 
for AA bilayer which can be derived at zero temperature 
from Eqn. (Ili 
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8 5(n) max(/i, 7) + 6[0 - 2\(i - j\] 

+ e[n- 2(^ + 7)]. 



(20) 



If the delta function here is replaced by a Lorentzian with 
broadening of 277, then this formula gives an excellent 
representation of the numerical results in Fig. [3] It is 
also possible to derive an expression for the imaginary 
part of the longitudinal conductivity which is Kramers- 
Kronig-related to Eqn. (|20|) by the relation 



a"(n) = -— P / 
* Jo 



o 2 -O 2 



(21) 



with ct'(uj), the real part of the conductivity given by 
Eqn. (HOI). Hence, 



00 
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= — max( 7 , fi) 
fi-2| 7 -/i| 



hi 



n - 2(7 + M ) 



f2 + 2(7 + /i) 



(22) 



If <5(f2) in the real part is replaced by a Lorentzian form 
as we have discussed, then instead of the Kramers-Kronig 
transformation of <5(f2) to l/(irQ,), we have r/[7r(f2 2 +r 2 )] 
transforming to Q,/[tt(VL 2 + T 2 )] where T = 2-q represents 
the transport scattering rate rather than the quasiparti- 
cle scattering rate 77 that enters the broadened spectral 
functions A(k,u). 

The issue of optical spectral weight redistribution with 
variation in chemical potential can be addressed globally 
by introducing the partial optical sum: 



duo. 



0+ CT o 



(23) 



which is defined as the area under the conductivity graph 
up to energy f2. For the real part of the longitudinal 
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FIG. 4. (Color online) Partial optical sum in units of 7 
for various values of chemical potential with (a) fi < 7 and 
(b) /j, > 7. Inset: The evolution of the positive frequency 
spectral weight W, found under half the delta function in the 
analytic solution for the conductivity, Eqn. (|20[) as a function 
of chemical potential. 



conductivity, is shown for various values of /1/7 in 
Fig. |U In all cases, at sufficiently high frequency the in- 
tegrated spectral weight returns to the /i = value (solid 
black curve of Fig. IUa)). The inset of Fig. @] shows the 
spectral weight of the delta function in the analytic solu- 
tion, Eqn. 1201 where we have taken only half the weight 
of the delta function as only half the function is present 
for positive frequency. The functions converge to the 
(i = case at frequencies above = 2(7 + /^) once most 
of the spectral weight from the Drude contribution is in- 
tegrated as well as the contribution from the interband 
edges. These curves once again provide an interesting dif- 
ferentiation between the /1 < 7 regime, where the Drude 
weight remains the same and the lower energy kink moves 
to lower J7, and the \x > 7 case where the low frequency 
part of the partial sum increases with \x and the low en- 
ergy kink moves to higher fi. In principle, such a quantity 
could allow for an experimental determination of 7 based 
on the transition from the behavior of one regime to the 
other. 

Turning to the perpendicular conductivity, we show 
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FIG. 5. (Color online) (a) Real part of the perpendicular 
conductivity for AA-stacked bilayer graphene for various val- 
ues of fx. The perpendicular conductivity is negligible every- 
where except near f2 = 2 7 where there is a sharp peak (note 
the axes of this plot). Here, the transport scattering rate is 
r = 2r/, where r\ = O.OI7 in all of our numerical work, (b) 
Real part of the perpendicular conductivity for AB-stacked 
bilayer graphene. A strong peak occurs 7 for finite [s, and an 
absorption edge occurs at 2 max(/i, 7). 



the response for the AA- and the AB-stacked bilayer in 
Fig. [5] The AA-stacked graphene has a strong absorp- 
tion associated with 27 which is finite at charge neutrality 
and increases with doping. It is also possible to derive a 
closed form algebraic formula for the perpendicular con- 
ductivity. For AA-stacking we obtain: 



c_lo 7 



(7-^) 2 ©(7-M)+27M 



(24) 



and by Kramers-Kronig transformation 



o±o 77r(f2 2 — 47 2 ) 



[( 7 -M) 2 e( 7 -M)+27M] 



(25) 



for the real and imaginary parts, respectively, where 
(7^o = (e 2 /4?i) • (7 d/Tivp) 2 . The physics of this case 
is as follows. For in-plane conductivity, charge carriers 
must hop from one sublattice to the other to produce a 



current, thus in the absorption process, interband tran- 
sitions are between two bands, each of which reflects the 
two sublattices by having different chirality label. Hence 
the transitions shown in the inset of Fig. [3fa) are be- 
tween two cones which have opposite chirality but the 
same bonding or antibonding wavefunctions. However, 
for the interlayer current in the AA-stacked case, the 
carriers hop between the A(B)-sublattice of one plane 
to the A(B)-sublattice of the other plane and as a result 
absorptive transitions for this form of transport will only 
occur between bands of the same chirality but different 
bonding which in reference to the inset of Fig. [3Ja) would 
be vertical arrows connecting the parallel bands in this 
case. As these arrows are always of length 27, there is 
only one very strong absorption peak at 27. 

Similar analytical results for AB-stacking at finite dop- 
ing are, to our knowledge, not in the literature and 
so we will provide it here for comparison (the equiva- 
lent form for AB-stacked in-plane conductivity has been 
given previously in the literature 3 -- 7 ). Some numerical 
work along with some analytical analysis has been done 
previously and our results are in agreement with those 
works. 36,59 In particular, Ando and Koshino considered 
polarization effects^ which we do not include here. The 
AB-stacked perpendicular response seen in Fig. [SJb) is 
quite different from the AA case. Absorption occurs at 
all frequencies but is on the scale of u±q. An absorption 
edge occurs at 27 similar to the AA case, although it 
is weaker by comparison and continues on to higher fre- 
quency as it results from transitions between the lowest 
and highest energy bands of the AB case in Fig. [5J which 
are the bonding and antibonding bands of the A-B dimer 
strongly coupled by hopping 7. 15 Moreover, absorption is 
seen at all frequencies and at finite doping a very strong 
peak occurs at 7 much as is seen in the in-plane con- 
ductivity for the AB-stacked bilayer—. It is not entirely 
surprising that the perpendicular conductivity displays 
elements of the in-plane conductivity with the same phys- 
ical origin^. The lower energy bands represent hopping 
between the non-dimer A and B sites in the two planes 
which must occur by first hopping in the plane to the 
neighbour site which is part of a dimer, hopping up the 
dimer bond and then over to the non-dimcrcd site in the 
second planed The AA-stacked case does not have this 
element. For AB-stacking, we derive an analytic formula 
for the perpendicular conductivity which is: 



<X ZZ (Q) 
C_L0 



where 











2(0 + 7) 2(fi- 7 ) 



9(0 - 2 7 ) 



9(0-2//) 



c( M )5(0- 7 ), (26) 



H(7 + H) 7, 2^ + 7 

c{fi) = -In 

7 2 7 



M(7 - M) , 7 ln 2/i-7 
7 2 7 



9(^-7). 
(27) 
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This formula also agrees quite well with the numerical 
work shown in Fig.[5jb) for various values of the chemical 
potential, provided the delta function is rewritten as a 
Lorentzian with broadening of T — 2rj. The expression 
for the imaginary part is given as: 



2il 



— 



1 



J2-7 
1 



-In 



2max(7, fj,) — Q, 



+ 



ft + 7 

27 



-In 



2fi + n 
2max(7, fi) + fi 



<>- 



-In 



r 



2^ + 7 



2max(/i, 7) — 7 
(28) 



where c(/i) is given by Eqn. (|27p. 

We can also examine the effect of adding a bias be- 
tween the two layers. To do this, we need to include an 
additional term in our Hamiltonian given by Eqn. (JT|) of 
the form^ 5 . 

H' =-A^2(a\ n a ln + b[ n+Si b ln+Sl ) 

n 

-2 A J2( a 2-n a ^ + bt n+sl b 2n+ sO- (29) 
n 

This bias raises the energy on the lower plane by +A/2 
and lowers the energy on the upper plane by —A/2 
providing an overall bias of A. For the case of AB- 
stacking, this introduces a gap in the energy dispersion 56 
and provides interesting features to the conductivity^ 5 -. 
For AA-stacked bilayer, this gives the energy dispersion 
e a (k) = ±[e + HV7 2 + A V4], where a = ±1. We 
can see that this is equivalent to a renormalization of the 
interlayer hopping parameter of the unbiased system to 
a value 7' such that 7' = -\/7 2 + A 2 /4 and will therefore 
introduce no new features into the conductivity. 



IV. CONDUCTIVITY OF AN AAA-STACKED 
TRILAYER 

These ideas can also be extended to trilayer graphene. 
For the case of AAA-stacked trilayer graphene, where A 
(B) sites in each layer arc stacked directly in line with the 
corresponding sites in the other layers, our Hamiltonian 
now becomes 



H 



in+S L 



■H.c. 



+ 7 E (4n°ln + b\ n bx n + 4„«2t. + &3„&2n + H.C. 
n 

(30) 

where i = 1, 2, 3 indexes each of the three layers. Here, 
we allow the usual nearest-neighbour intralayer hopping 
t and interlayer hopping 7 between the neighbouring 
planes; again, we have ignored hopping from an A(B) 



site in one layer to a B(A) site in another layer as well as 
hopping from an A1(B1) site to an A3(B3) site as these 
hopping energies are very small 5 - 5 ". Transforming to k 
space, we obtain the following matrix representation: 



H 



( 











7 













7 


f*(k) 


7 











/(fc) 


7 








7 


f (fc) 













f(k) 


7 











V /*(*) 


7 











/ 



(31) 



where we have used the eigenvector ^ = 
(aik,b 2 k,a3k,b 3 k,a2k,bi k ). Reflecting the fact 
that we now have six atoms per unit cell, we obtain the 
following six energy bands: 

s(k) = ±|/(fc)|,±[|/(fc)| - v^y],±[|/(fc)| + V27], (32) 

where |/(fe)| is the energy dispersion of monolayer 
graphene, equal to e = Tivpk at low energy. The first two 
bands, ±e, are the original graphene bands, the second 
two bands, indexed ±£i(fe), and final two bands, indexed 
±e-2{k), are monolayer bands shifted by =Fv / 27, respec- 
tively. A plot of the band structure can be seen in the 
inset of Fig. [51 From this, we see that the trilayer is like 
the sum of a monolayer and a bilayer with an interlayer 
hopping of y2~7- 

Using the same formalism as before, we can derive an 
expression for the real part of the zero temperature lon- 
gitudinal conductivity in terms of the spectral functions; 
we obtain: 



2n 



fi—Q 



2^ 



d 2 k 

(27T) 



2 2v% 



x [2A n (u + Q)A n (u>) + 2A 13 (uj + Q)A 13 (l>) 

+ 4A 15 (uj + n)A 15 (uj) + A 22 {uj + n)A 22 (uj)}, 

(33) 



where 



A U = - 



A 13 = T 



2S(u-e) + 26(u + e) + 5(uj - £l ) 

+ S(U3 + £1) + 5(u> - £ 2 ) + S(uj + £ 2 ) , 

(34) 

- 25(u-e) - 26(lu + e) + S{uj - £1) 

+ 5(lu + £1) + 5{uj - £ 2 ) + 5(uj + £ 2 ) 

(35) 



A 15 = 



V2ir 



S(uj - £1) + S(lo + £1) + 5(uj - £ 2 ) 



5(uj + £ 2 ) 
(36) 
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FIG. 6. (Color online) Conductivity for AAA-stacked tri- 
layer graphene. Interband absorption edges appear at SI = 
2/Li, 2|%/2~7 — /x| and2(%/2~7 + /i). These transitions are shown 
on the band structure given in the inset. 



and 



A 22 = - 



S(lu - e\) + 5(u + ei) + 5(uj - e 2 ) 



+ 5(lu + e 2 ) 



(37) 

Several numerical curves for the conductivity can be 
seen in Fig. [6l For finite //, there are three steps in the 
conductivity each of value <To, leading to a constant back- 
ground at high frequency equal to three times that of 
a single sheet, reflecting the trilayer nature of the sys- 
tem. In each case, there is always an absorption edges at 
fi = 2/i. The other two edges occur at 2\/i — y/2j\ and 
2(/i + v27), where the former decreases with increasing 
fi for fx < y^2j and then increases for /i > V2"f, while 
the latter always increases with /i. This is a similar pat- 
tern to the AA-stacked case and so the combination of 
the above confirms that the trilayer acts as the sum of 
a monolayer plus bilayer. Even at fi = 0, the flat back- 
ground of (7 , due to monolayer behavior, is added to the 
unusual Drude plus 2cto interband behavior of the fi = 
bilayer shown in Fig. [31 but for effective interlayer hop- 
ping of V27- The ability to tune the flat background in 
the IR spectral region from — > 3cto in steps of oq by 
changing the doping is an interesting feature that could 
possibly be of some advantage to technological applica- 
tions. 

Given Eqn. (1331) . an analytical expression for a xx of 
AAA-stacked trilayer graphene may be written down. It 
has the expected form: 



c(fi) 



(TO 



=4<5(fi)[2max(/z, ^2 7 ) + fj] + 0[fi - 2/x] 

+ 6[fi - 2(/i + V2 7 )] + 6[fi - 2|/i - V2 7 |], 

(38) 



which stresses the existence of three decoupled and 
shifted monolayer dispersions where interband transi- 
tions are only permitted between the corresponding cones 
shown in the inset of Fig. [51 

The imaginary part is again found by applying 
Eqn. ([211) to Eqn. (08]). It is 

SI - 2/i 



= 4-[2max(/i, V27) + fi] + -In 

(To Mi 7T 



7T 



fi - 2(/x + V2ry) 



fi + 2(/i + V27) 



Im 

7T 



fi + 2/i 

fi - 2(/i - y/2~f) 



fi + 2(/i - ^2 7 ) 



(39) 



V. AA-STACKING WITH SPIN ORBIT 
COUPLING 

Finally, to cover a variety of possible scenarios, we ex- 
amine the effect of spin orbit coupling (SOC) on AA- 
stacked bilayer graphene. Such effects in AA- and AB- 
stacked bilayers were considered by Prada et in the 
context of studying systems which may manifest a topo- 
logical insulating phase. These authors have studied both 
the case of SOC in each plane and SOC in only one plane, 
the latter case taken to be a toy model for spin-orbit 
proximity effect. For our purpose here, these consider- 
ations illustrate the generic features of opening energy 
gaps in the AA-stacked band structure at k = and at 
the charge neutrality point. Recall from our previous dis- 
cussion that biasing the bilayer does not open a gap in 
the AA-stacked case as it does in the AB-stacked case, 
but SOC will do so. 

For a single sheet of graphene including SOC, the tight 
binding Hamiltonian is 60 



At z s z f(k) 
f*(k) -At z s z 



(40) 



where, in the continuum limit, f(k) = h>VF(r z k x — iky") 
and A = 3^/3t so , with t so , the next-nearest neighbour 
hopping amplitude. t z — ±1 for the Dirac points K 
and K' and s z = ±1 corresponding to the up/down spin 
component perpendicular to the graphene sheet*^. With 
this, the Hamiltonian for the AA-stacked bilayer is 



n lT„f 



(41) 



where we have used the eigenvector \1/ = 
(«ifc, b\k, a 2 k, b 2 k)- When dealing with the case of 
SOC in both layers, both I/j 1 ^ ^ and i^r^s* are given 
by Eqn. ([4H| and H± is the coupling between the layers, 
again taken to be 

*X = (J°) (42) 
for AA-stacking. Eqn. (|41[) gives the four energy bands 43 



e a (k) = ± Ve 2 + A 2 + (- 



(43) 
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where a = 1 and 2. These are illustrated in the inset of 
Fig. (7{b) for a particular K point. 
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FIG. 7. (Color online) Conductivity Om^) versus fi/7 for 
the case with SOC in both layers, (a) fi = and A/7 = 0, 
0.75 and 1.2 in each layer, (b) The case of finite doping with 
(i = O.57 and A = 0.757 with that for A = shown for 
comparison. The inset shows the band structure for this case 
and the transitions that occur to give the absorption edges 
seen in the main frame. 

While we are interested in a xx , we will forgo provid- 
ing the explicit details of the calculation as they can be 
developed by following the procedure already outlined 
earlier. The expression for a xx (fl) given in Eqn. (|10j) can 
still be used. However, in lieu of the degeneracy factor 
Nf, a sum over t z and s z should be taken when using 
Eqn. (|4ip to calculate the Green's function. Our velocity 
operator v x written in the basis used for this section can 
be evaluated as before from Tiv x = dH/dk x . 

The effect on the real part of the longitudinal con- 
ductivity of AA-stacked bilayer graphene when SOC is 
present in both layers can be seen in Fig. [7] for the case 
of (a) /1 = and (b) /1 = O.57. The two shifted monolayer 
dispersions in the band structure are now gapped by 2A 
about ±7 (see inset of Fig. [7])- For fi = and A < 7, as 
shown in Fig. EJa) , there is the usual Drudc conductivity 
and a jump at 27 as in the case of no SOC. However, the 
shape of this jump is typical of a gapped electronic spec- 



trum which gives rise to a discontinuity in the electronic 
density of states. Indeed for the monolayer of graphene, 
such behavior has been calculated.— As the frequency 
increases the usual bilayer background is recovered. The 
features of the conductivity show that, as in the case of no 
SOC, transitions are only allowed within each decoupled 
(like-colored) monolayer band. For the case of A > 7, 
we no longer have the Drude contribution as no states, 
and therefore no transitions, are available at zero energy. 
The peak in the conductivity now occurs at 2A and again 
we can see that we obtain the usual background conduc- 
tivity for significantly high frequency. We have chosen 
our values for A/7 to be in keeping with the parameters 
used by Prada et alM-, however, in graphene the intrin- 
sic SOC gap is ~ 10~ 3 meV and hence A/7 ~ 10~ 5 . 
Nonetheless, we have chosen this model to indicate the 
effect of energy gaps appearing in the band structure. 
Indeed, other graphene-like systems are now being stud- 
ied which have much larger predicted SOC gaps, such 
as silicene (A ~ 1.5 meV) and germanene (~ 25 meV) 
which can be further tuned by strain or perpendicular 
electric field . 61 ' 62 Likewise, a varying mass gap of up to 
~ 150 meV has been found in the 3D topological insula- 
tor TlBi(Si_ x Se x ) 2 by Sato et al. 63 and the SOC gap in 
monolayer M0S2 and other group- VI dichalcogenides is 
on the order of 1.5— 1.8 eV— . Consequently, the results of 
this section may be very relevant to future developments 
in these graphene-like systems. 

In Fig. [7jb), the effect of finite doping is considered in 
comparison with the case of no SOC. For small enough 
A, the SOC curve will track the one without SOC with 
the exception that there is a peak at each absorption edge 
reflecting the energy gaps in the band structure at k = 0. 
If A > I7 — fi\, the edge will be at 2 A rather than at 
2I7 — /1 1, and likewise for A > 7+/X, there will be only one 
jump at 2A. A Drude contribution remains provided that 
A < 7 + fi. The behavior of the absorption with doping 
reflects possible transitions between the original shifted 
monolayer bands subject to the opening of a gap. Mixing 
between these two monolayer-type bands does not occur. 
This is illustrated in the inset of Fig. [TJb) . The behavior 
embodied by these figures can be derived analytically and 
we find it to be: 



a xx (Q) = a xx (fl, I7 - fi\) + 5-^(0,7 + fi), 

with 

a xx (n,T) „(T 2 -A 2 ) 



(44) 



on 



T 



-<5(fi)9(T-A) 



2A 
7T 



0[ft-2max(T,A)] 



(45) 



where this last expression is the conductivity for mas- 
sive Dirac quasiparticle o 26 ' 65 . This formula is in good 
agreement with the numerical work. Note that for the 
longitudinal optical conductivity, there is a conservation 
of spectral weight upon introducing a finite [i and finite 
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A. The corresponding imaginary conductivity is given 
by 

= 5^(n, | 7 - Ml) + 5^(«.7 + m), (46) 

where 



4 T 



ttO T 
4A 2 r 

+ 



-e(T 
1 



A) 



1. 



-In 



Q - 2max(T, A) 



1 



max(T, A) SI 



-hi 



Q + 2max(T, A) 
- 2max(T, A) 



2max(T,A) 
(47) 



If we only keep SOC in one layer, one of the diagonal 
elements of our Hamiltonian given by Eqn. (|41[) becomes 



f(k) 
f*(k) 



(48) 



and the four energy bands becomi 



,43 



£ a {k) 




A 2 

+ 7 2 + ^r + (-) Q 



A 4 

4e 2 7 2 +7 2 A 2 + — , 



(49) 

where a = 1 or 2. The conductivity for this case can be 
seen in Figs. |8] and |9] for the case of fi = and finite n, 
respectively. The band structure is plotted in the insets 
of Fig. [SJa) and (b) for A = O.57 and A = I.27, respec- 
tively. The key energy levels labelled in the insets are 
given by 



\ 



?2 



A 2 . 
- + ( - 



'4e 2 + A 2 



A 4 

T' (50) 




FIG. 8. (Color online) Conductivity for the case of SOC in 
one layer only and the doping set to charge neutrality. The 
value of A in the one layer is (a) A = O.57 and (b) A = I.27. 
The band structure for each case is shown as an inset. 



where A = A/7 and e = e/7. Gaps now appear about 
zero energy and about ±7. 

The shape of the optical curves reveals that transitions 
are now allowed between every band unlike in the pre- 
vious case with SOC in both layers. In Fig. a strong 
absorption is seen at 2E\ 1 (see label in inset to the figure) 
which is ~ A at small values of A. The sharpness and 
strength of this absorption feature is due to the existence 
of a square root singularity in the electronic density of 
states at half this energy or ~ A/2. This feature has 
replaced the Drude absorption due to the gap about zero 
energy, but very little else changes in Fig. [8fa) from the 
no-SOC case for small A < 7. For A > 7 shown in 
Fig. E[b) , further structure appears which can be traced 
to various transitions in the band structure as has been 
labelled in the inset. While the behavior is relatively 
simple for \i = 0, for finite fi it is much more compli- 
cated. Fig. [^a) and (b) shows the case for A = O.57 
and A = I.27, respectively, where fj, is varied through 
key parts of the band structure [refer to the insets of 
Fig. [5] . While the high frequency behavior is similar to 
what we have seen for no SOC with an interband absorp- 
tion edge tracking 2(7+//), the low frequency behavior is 



very structured. Notable is that the feature at 2Ei 1 (or 
~ A for small A) is quickly suppressed by Pauli block- 
ing at finite /i and there appears to be a deep absorption 
minimum just on or after 17 = 27. Comparing this figure 
with Fig. [3] shows that for each value of /j, the underly- 
ing structure of the double step for the AA-stacked case 
with A = is retained, while the low frequency behavior 
oscillates about the A = case. The governing behavior 
still appears to be dominated by the decoupled mono- 
layer bands however the complicated structure arising 
from various interband transitions within the different 
bands is reminiscent of what is found in Bernal-stacked 
bilayer graphene for finite \x and an asymmetry gap2£ 
which gives rise to similar looking (but not quite the 
same) band structure. We do not have an analytic form 
for this complex behavior and one must rely on an ex- 
amination of the optical transitions available in the band 
structure to identify the detailed structure. However, it 
is clear that the case of SOC in only one plane is quite 
different from that where SOC is in both planes. 
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FIG. 9. (Color online) As for Fig. [8] but for finite doping as 
indicated in the legends. The SOC is only in one of the two 
layers with value (a) A = 0.57 an d (b) A = 1.2f. 



an important energy scale. In contrast to the monolayer 
which is constant at all frequencies for charge neutrality, 
the bilayer exhibits a Drude conductivity and absorption 
is Pauli-blocked for frequencies less than 27. At finite 
doping, a double step occurs in the interband absorption 
with onset for each step at I7 — fi\ and 7 + ji and the 
behavior is non-monotonic in fi as a result. The perpen- 
dicular response is completely centered on 27 with the 
absorption peak being very significant at all dopings and 
increasing with doping. This is unlike the behavior seen 
in AB-stacked bilayer graphene. Applying a bias across 
the bilayer does not open an energy gap in the band struc- 
ture but merely renormalizes the effective interlayer hop- 
ping to a greater value. The conductivity of the trilayer 
exhibits three absorption edges leading to a flat conduc- 
tivity background which steps up through <7o to 3<tq . The 
conductivity in this case is clearly seen to be a sum of 
that due to a monolayer added to that of a bilayer with 
interlayer hopping \/2j. Including spin orbit coupling in 
each layer of the AA-stacked bilayer leads to double step 
behavior where the absorption edges show a peak due to 
a gap of 2A in the band structure and the location is set 
by an interplay between the energy scales of chemical po- 
tential [a, 7 and A. With spin orbit coupling in only one 
plane, an absorption peak is found at about A and very 
complicated structure as a function of frequency results 
at finite /1 although there remains a remnant of the un- 
derlying structure of simple AA-stacked bilayer conduc- 
tivity. As experimental isolation of AA-stacked graphene 
has been reported and interest in topological insulators 
with spin orbit coupling is high, this study may be timely 
for future work in these areas. 



VI. CONCLUSIONS 
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